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In order to examine the gravitational waves emitted Irom the neutron stars in the tensor-vector- 
scalar (TeVeS) theory, we derive the perturbation equations for relativistic stars, where for simplicity 
we omit the perturbations of vector field. That is, we consider the perturbations of scalar and tensor 
fields. With this assumption, we find that the axial gravitational waves, which are corresponding to 
the oscillations of spacetime (lo modes), are independent from the perturbations of scalar field and 
the effects of scalar field can be mounted only via the background properties. Using two difi'erent 
equations of state, we calculate the complex eigenfrequencies of axial w modes and find that the 
dependences of frequencies on the stellar compactness are almost independent from the adopted 
, equation of state and the parameter in TeVeS. Additionally, these dependences of frequencies of 

' axial u) modes in TeVeS is obviously different from those expected in the general relativity. Thus 

\ the direct observations of gravitational waves could reveal the gravitational theory in the strong-field 

regime. 
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I. INTRODUCTION 

I , As an alternative gravitational theory, the tensor-vector-scalar (TeVeS) theory has attracted considerable attention. 
^jPj' This theory is proposed by Bekenstein [l[ as a relativistic theory of Modified Newtonian Dynamics (MOND) [2], and 
TeVeS reproduces MOND in the weak acceleration limit. The most important advantage to adopt TeVeS might be 
the point to be possible to explain many galactic and cosmological observations without the need for dark matter 
Q . In TeVeS, one has successfully explained the galaxy rotation curves and TuUy-Fisher law without the existence of 
^ '^ 1 dark matter T|. Additionally, TeVeS is possible to explain not only the strong gravitational Icnsing but also the 
I ' galaxy distribution through an evolving Universe i(| without cold dark matter. On the other hand, in the strong-field 
T^lj- . regime of TeVeS, Giannios found the black hole solution by solving the field equations for static, spherically symmetric 
C^l ' spacetime in vacuum i.e., this corresponds to the Schwarzschild solution in the general relativity (GR), where he 
, found two distinct branches of solutions dependent on the form of the vector field. Subsequently, Sagi and Bekenstein 
' generalized the Schwarzschild solution in TeVeS to the Reissner-Nordstrom solution 0], while Lasky, Sotani, and 
^\ ' Giannios derived the Tolman-Oppenheimer-Volkoff (TOV) equations and they produced the neutron star models in 
O TeVeS '^]. 

^ ] Recently, it was discussed the possibilities of distinguishing TeVeS from GR. For example, it was suggested that 
• ^ . one can distinguish the gravitational theory with the redshift of the atomic spectral lines emanating from the surface 
^\ ' of neutron stars or with the observations of neutron star oscillations via the emitted gravitational waves p^ . 



Additionally, Desai, Kahya, and Woodard pointed out another test of TeVeS where they showed that there is an 
appreciable difference in the Shapiro delays of gravitational waves and photons or neutrinos from the same source and 
suggested the possibility to test of TeVeS by observing this difference emitted from gamma ray bursts or core-collapse 
supernovae. In this article, we examine whether observations of gravitational waves emitted from the neutron stars 
can provide an alternative way of probing the gravitational theory in the strong-field regime, where we focus especially 
on the gravitational waves associated with the oscillation of spacetime itself. 

In fact, the tests of gravitational theories in the strong- field regime are quite important, because the gravitational 
theories in the strong-field regime are still largely unconstrained by observations in contrast to those in the weak-field 
regime, which have been subject to numerous experimental tests. However, owing to the developments in observational 
technology, it is becoming possible to observe compact objects with high accuracy. Via the observations of X-rays 
and/or gamma rays emitted from compact objects, one can know the properties of compact objects and could use as a 
direct test of the gravitational theory in the strong-field regime. As an alternative way to observe compact objects, the 
gravitational waves are also expected to obtain the raw information of the compact objects. With these observational 
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properties, it might be possible to distinguish the gravitational theory in the strong-field regime [13. [T4|. 

Meanwhile, it has been suggested that the observations of gravitational waves can provide a unique tool not 
only for estimatin g th e stellar parameters such as mass, radius, rotation rate, magnetic field, and equation of state 
(e.g., [l^ [l^, [TtI Tisl [T9I I^O*!), but also for verifying the gravitational theory, which is called "gravitational- wave 
asteroseismology" . Furthermore, the detailed analysis of the gravitational waves also makes it possible to determine 
the radius of accretion disk around supcrmassive black hole [2l|] or to know the magnetic effects during the stellar 
collapse [22!|. The eigenmodes, which are mainly excited during the formation of a neutron star or during the starquakes 
and emit detectable gravitational waves with the ground-based gravitational wave detectors such as LIGO, GEO600, 
VIRGO, and TAMA300, are the fluid / and p modes and the w modes [l^l, which are associated with oscillations of 
the spacetime [2J| . The possibility to distinguish TeVeS from GR by using the fluid oscillations (/ and p modes) has 
been already discussed in [To| . where the Cowling approximation was adopted. Thus in this article we focus on the 
w modes. The w modes are similar to quasinormal modes of black holes. They have higher frequencies and shorter 
damping times than the fluid modes, i.e., in GR the typical frequencies are around 7-12 kHz and damping times 
are order of 0.1 ms. In general, the oscillations on the spherically symmetric spacetime can be classified as axial and 
polar depending on their parity. Since the axial w modes are known to have the same qualitative behavior as the 
polar w modes in GR, in this article we will examine only the axial w modes. Moreover, for simplicity, we consider 
the perturbations of scalar and tensor fields while those of vector field will be omitted in this article. The more 
complicated analysis of the polar w modes and/or including the perturbations of vector field will be seen somewhere. 

This article is organized as follows. In the next section, we describe our notation and briefly introduce the theoretical 
framework of TeVeS. In Sec. Illll we derive the perturbation equations for the axial perturbations. Then the oscillation 
spectra of neutron stars in TeVeS are shown in Sec. llVi finally we discuss the results related to gravitational wave 
asteroseismology in Sec. |Vl In this article, we adopt the unit of c = G = 1, where c and G denote the speed of light 
and the gravitational constant, respectively, and the metric signature is (— , +, +, +). 



II. STELLAR MODELS IN TEVES 



In this section, we mention the stellar models in TeVeS. TeVeS is based on three dynamical gravitational fields, 
such as an Einstein metric f^^^, a timelike 4- vector field lA^ , and a scalar field (p, in addition to a nondynamical scalar 
field (7. The vector field fulfills the normalization condition g^i^U^W^ = —I and the physical metric is given by 

5m- = e""'^ (5m- + ^M^-) - e^'^U^U,. (2.1) 

All quantities in the physical frame are denoted with a tilde, and any quantity without a tilde is in the Einstein frame. 
Varying the total action with respect to g^"^, one can obtain the field equations for the tensor field 

+ e^„ (2.2) 

where is the energy-momentum tensor in the physical frame, Ta(^^Ui,^ = Ta^Un + Ta^U^, and G^^ is the Einstein 
tensor in the Einstein frame, while the other sources r^i/ and O^i^ are given by Eqs. (2.4) and (2.5) in [l3|. Notice 
that the conservation of energy-momentum is given in the physical frame as V^T'^'' = 0. Additionally, by varying the 
total action with respect to and (p, one obtains the field equations for the vector and scalar fields as Eqs. (2.6) 
and (2.7) in [l^ . which include two positive dimcnsionless parameters, k and /C. These are the coupling parameters 
for the scalar and vector fields respectively. 

A static, spherically symmetric metric in Einstein frame can be expressed as 

ds^ = g^pdx^'dx'^ = -e'^^^'^dt^ + e'^^^'^dr'^ + r^d^"^, (2.3) 

where c?ri^ — dO^ + s\v?' BdcjP' and = 1 — 2TO(r)/r, while the vector field can be given as Z//^ — {U^{r),U^{r),0,0). In 
this article we set W=0, because it has shown that in vacuum, the parameterized post-Newtonian (PPN) coefficients 
for a spherically symmetric, static spacetime with a non-zero can violate observational restrictions Q. Then with 
the normalization condition one can show that = (e"''/^, 0, 0, 0). With this vector field, the physical metric is 

ds^ = go.pdx^dx'^ = -e^+^fdi"^ + e'^-^'^dr^ + e-^'^r^dfl^. (2.4) 
Finally we assume that the stellar matter consists of a perfect fluid 



Gai/ = SttG 



T„. + (l-e-4^)Z^"f„(^W,) 



Tfj^u = [p + P\ Uf,u^ + Pg^^, 



(2.5) 
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where m^, p, and P are the four- velocity of the fluid, the total energy density, and the pressure in the physical frame. 
The concrete neutron star models in TeVeS can be seen in the previous papers [1, • 

To construct the stellar models in TeVeS, three new parameters, fc, K and (/5c, are introduced with respect to GR, 
where ipc denotes the cosmological value of the scalar field. The value of k has a tightly constraint as fc 0.03 by 
both cosmological models and also planetary motions in the solar system, while the cosmological considerations imply 
that the value of Lpc is restricted to < (/Sc <C 1 With respect to the value of (pc, it was shown that ipc could have 
a minimum value of around 0.001, which is based on the causality issues inside the star Anyway, since it was also 
found in 9] that the neutron star models are almost independent from the values of k and ^pc^ in this article we make 
examinations with k = 0.03 and Lpc = 0.003. On the other hand, although we pointed out in the possibility to 
determine the value of IC with the observations of gravitational waves related to the stellar oscillations, so far there 
is no severe restriction on /C except that K, should be in the range of < /C < 2 [1, [§] . Thus in this article we adopt 
various values of K, in the range of < /C < 2 and study the dependence of gravitational waves on the parameter K,. 
Furthermore, as equilibrium stellar models, in this article we adopt the similar models in i.e., two different EOSs 
such as EOS A and EOS II. 

At the end of this section, we introduce two useful properties, i.e., the total Arnowitt-Deser-Misner (ADM) mass 
-^^ADM and the scalar mass , which are defined as 

Madm = fm, + ^^^'j (2-6) 



47r 

M^ = A7T I r^(p + 3P) e'^'^+'^^/^-^^dr, (2.7) 



^ ' 

where rric is the mass function evaluated at radial infinity. Notice that the scalar mass is constant outside the star, 
while with ADM mass one can describe the asymptotic behavior of the physical metric as 

+ ^ + (2.8) 
^„ = l + ^^ + ofl), (2.9) 

where f = re~'^'' . 



III. PERTURBATION EQUATIONS 



As mentioned before, for simplicity we assume that 5Ui = in this article, where the subscript i corresponds to 
i = r, 9, and </>. In fact there is no physical reason for this assumption. But in order to construct the neutron star 
models with the observed masses, the value of the coupling constant for vector field, IC, should be smaller than at 
most IC ~ 0.8, which depends on the adopted EOS M- So since one could expect that with smaller value of IC 
the effect of the vector-field perturbation on the spacetime oscillations might be small, in this paper as a first step we 
neglect this type of perturbation to make our problem simple. (Also see in Appendix [X] for the effect of the vector-field 
perturbation on the spacetime oscillations.) It is noteworthy that SW^ ^ and 5Ut ^ even with the assumption 
that SUi = 0, which are results from the metric perturbations. Under this situation, the perturbations of the scalar 
and vector fields are expressed as 



Sip = 6ip{t, r)Yirn, 

'1 



<^5tt, 0,0,0 



while, using the Rcggc- Wheeler gauge, the perturbed metric tensor in the physical frame is given as 



(3.1) 
(3.2) 

(3.3) 



where denotes the axial part of metric perturbations 
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and Ti'i^i) denotes the polar part of metric perturbations 
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r'^Kim sin^ i 



(3.5) 



Here, the functions /io,;m, ^i,im, Ho^im, i?i,imj ^^2,;™, and ifjm describing the spacctime perturbations have only 
radial and temporal dependence while Y/m = Yi^id, 4>) is the spherical harmonic function. Since, from Eq. (|2.ip . the 
perturbed metric tensor in the physical frame is connected to that in the Einstein frame 6g^,^ as 

= -2e-2^ (5^^ + Uf,U^) Sip + + + Z^^^Z^^) - 2e^m^JJi,Sf - e^'^ {5U^,W + U^,5U^) , (3.6) 

the perturbed metric tensor in the Einstein frame is 

/ {e-^'^Ho.im + 26^) e^^H^.i^ -e^'^/io,;™ sin"! 68^ 

^ -e^'^/ii,;™ sin-i 05^ 



= E E 
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2Sip) sin^ ( 



Y, 



Im ■ 



J 



(3.7) 

This expression of Sgf^i, is used when the perturbations in the Einstein frame will be transformed back to the physical 
frame. Now, by defining the new set of perturbation functions, Ho,im, -f^i.Zm, ^^2,;™, Ki„i, /lo.im, and hi^im, as follows 

Ho,tm = e-^-'^Ho^im + 26ip, (3.8) 

e^^Hi^i^, (3.9) 

e^^H2j„r + 25ip, (3.10) 

e-Ki„, + 25^, (3.11) 

e^^ho.im, (3.12) 



H2Um 
Khn 
ho,lm 



Im 1 



(3.13) 



the perturbed metric in the Einstein frame is simplified considerably and reduced to the "standard" Regge- Wheeler 
form of perturbed spherical metric. Here we emphasize that the scalar perturbation 5(p is linked only with the polar 
perturbations Ho^im, i?i,;m, H2.im, and Kim, while the axial perturbations ft-o.im and hijm are affected only by the 
contribution of the scalar field to the background. 

The perturbation equations for the gravitational waves can be obtained by taking the variation of tensor field 
equations (|2.2p . In order to derive the perturbation equations, we define the variations of pressure and energy density 
in physical frame as 



SP = SPYlm, 

Sp = SpYim, 

while the variation of four-velocity in physical frame as 

1 



6u' 
Sit''' 



-Zv-y/2jjY, 



1 



~^~vl2 



WYi„ 



-'^-"/^ VdgY 



smf 



^dsYir 



sin^ I 



-ip-ii/2 



{Vd^Yun + us\nedeYi„i) 



(3.14) 
(3.15) 

(3.16) 
(3.17) 
(3.18) 
(3.19) 



where (5P, 5p, W, V, and u are functions of t and r. Using these definitions, from the tO, t(j)^ r9, and r(j) components 
of linearized Einstein equations, we get 



1 ™ ^ 



(J = 0,1), 



l,m 



Im 



smt 



d^Ylra \=Q ( J = 0, 1) 



(3.20) 
(3.21) 
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where the explicit expressions of and f3\'^ are given in Appendix |Bl Notice that are some hnear combinations 
of polar perturbation functions, while (3\^ consist of only axial perturbation functions. Furthermore, from the 9(j) 
component and the subtraction of 09 and (pcf) components, one obtains two more equations 

^ {simXlm - thn aul9Wlm} = 0, (3.22) 

l,m 

J2 {tl^^^hn + Shn sill9Wlm} = 0, (3.23) 
l,m 

where sim and tim describe polar and axial type perturbations, respectively (see Appendix [B| . while Xim and Wim 
are defined as 

/ cos 9 \ 

Xi„,^2d^idg~—-)Yi,n, (3.24) 
\ smt^ / 

Wi„. ^ (dl - '-^do - -A-A) Yi^. (3.25) 
\ smt^ sm^ 

Additionally, from the tt, tr, rr components and the sum of the 99 and <j)(l) components, one can get four more 
equations 

^A™^i™=0 (7 = 0,1,2,3), (3.26) 

but we do not care these expressions in this article since these equations are corresponding to polar perturbations. 

Then, by taking the product of Eqs. p.20p - p.23p and (I3.26|) with Yim, integrating over the solid angle, and 
paying attention to the fixed values of I and m, we get ten partial differential equations in the variables t and r as 

A™=0, «L?=0, m.n = Q, (3.27) 
A™=0, i,,„ = 0, (3.28) 

where / = 0, 1,2,3 and J = 0, 1. Here Eqs. p.27p describe the polar perturbations, while Eqs. (j3.28|) describe the 
axial perturbations. It is worth noticing that the analytic expressions for Eqs. (|3.28p . i.e., Eqs. (|B3p . (|B4p and (|B6p . 
do not involve the perturbation of scalar field Skp. Thus the scalar perturbation is coupled only to the gravitational 
waves with polar parity. 

Combining Eqs. (|3.28p . one can easily derive a wave equation for the axial perturbations as 

$ - e^"-'^)/^ (^e'^'^-0/^^''^\ e"-^ y" + ^ l^j,' _ ^' + _ JLgC 1^^' + _ + i)^ + Z^i^GPr^ e'^^'^ $ = 0, 

(3.29) 

where the new function <i>(t, r) is introduced, which is defined as h\ = e('»^'')/^r$. Note that if /c = /C = this 
equation reduces to the standard wave equation describing axial perturbations in GR (25l |. As mentioned the above, 
the wave equation p.29p does not include the perturbations of scalar field and the effects of scalar field will enter 
only via the background properties. That is, the axial gravitational waves can be studied independently neither from 
polar gravitational waves nor from the scalar field perturbations. This axial gravitational waves are well-known as w 
modes (23|, which are quasinormal modes describing the pure spacetime oscillations and similar to the quasinormal 
modes of black holes. 



IV. SPACETIME PERTURBATIONS IN TEVES 



In order to determine the quasinormal frequencies of the axial w modes, one can think up two different techniques. 
The first approach is the direct time evolutions of Eq. (|3.29p and the frequencies will be determined by computing 
the Fourier transform of the signal at infinity, while in the second approach, one assumes a harmonic time dependence 
of the perturbations and the frequencies will be obtained by solving the eigenvalue problem with the appropriate 
boundary conditions. The first approach is quite simple, but with this approach one can identify only those of the 
quasinormal modes which are excited significantly, and the outcome depends strongly on the choice of the initial data. 
Additionally, using time evolutions it is quite difficult to identify quasinormal modes that damp out very fast. On 
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the other hand, although the second approach is more comphcated than the first one, it is possible to calculate both 
slowly and strongly damped quasinormal modes. Thus, in this article we adopt the second approach to determine 
the quasinormal frequencies of axial w modes. Assuming a harmonic time dependence as r) = <I>(r)e*"*, the wave 
equation (|3.29p can be rewritten as 



2 C-i/ 



— I' 
2 



f + 5 

r 



-^e'^ |rm' + 7m - l{l + l)r + 32TrGPr^e-^'^^ 



$ = 0. (4.1) 



Then, with appropriate boundary conditions, the problem to solve becomes an eigenvalue one with respect to the 
complex quasinormal frequencies lu. Physically, the real and imaginary parts of the complex frequency are corre- 
sponding to the oscillation frequency and the damping rate of each eigenmode, respectively. The imposed boundary 
conditions are that $ should be regular at the stellar center and that there are no incoming waves at infinity. In 
fact, near the stellar center, one can show that $ has a behavior of the form as <l>(r) = <i>or'+-'^ (l + 0{r^)), where 
$0 is some arbitrary constant. With this boundary condition, we can just integrate the above differential equation 
inside the star. On the other hand, outside the star, we use appropriate asymptotic expansions to meet the boundary 
condition at infinity, where we adopt the Leaver's continued faction method 26]. The concrete numerical procedure 
is described in Appendix [Cl 

In Fig. [l] as examples, we show the eigenfrequencies of w and wu modes for neutron star models with Madm — 
IAMq, where the open symbols are corresponding to the wu modes and the solid ones denote the w modes. From 
the observational point of view, the lowest w modes might to be relevant for the gravitational wave detectors [23| , 
while the higher w modes, i.e., W2, ws, W4, ■ • ■ , are probably difficult to detect if not impossible. As regarding to 
the Wu mode, the oscillation frequency is also suitable for detection by ground-based interferometers, although they 
damp out very fast. 
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FIG. 1: The complex frequencies of the axial w and wn modes for / = 2, where the stellar masses are fixed to be Madm = IAMq . 
The left panel corresponds to EOS A and the right panel to EOS II. In the figures, wu modes are shown with the open symbols 
and w modes with the solid ones. The circles denote the frequencies in GR, while the other symbols denote those in TeVeS 
with different values of /C. 

In the previous studies in GR, it is well known that the frequencies of w and wu modes can be described as a 
function of stellar compactness Madm/R p^[23l|. Similarly, we will examine the dependence of frequencies of wu and 
w modes on the stellar compactness. Figs. [5] and [3] show the frequencies of axial wu and wi as functions of the total 
stellar compactness, where the solid symbols are results with EOS A and the open ones are those with EOS II. From 
Fig. [21 we can see that the dependence of frequencies of wu modes on the stellar compactness are almost independent 
from the values of /C and the adopted EOS. Additionally, it is found that the deviation from the frequencies expected 
in GR is very little. We can see little difference between the expectations in GR and in TeVeS in the damping rates 
(imaginary part of complex frequencies) for the stars with larger compactness and in the oscillation frequencies (real 
part of complex frequencies) for the stars with weaker compactness. For the oscillation frequencies in TeVeS, we can 
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get the empirical formula, such as 



Re (wMadm) = -0.0412 + 0.7084 



M 



•'-'ADM 



3.305 



ADM 



R J \ R 

and the oscillation frequencies can be in very good agreement with this empirical formula. 



(4.2) 
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FIG. 2: For wu modes, the damping rates (left panel) and the oscillation frequencies (right panel) as functions of stellar 
compactness Maom/R- The solid symbols correspond to the results for EOS A and the open ones to those for EOS II. 



On the other hand, in Fig. [3] for the frequencies of wi modes, we can observe same feature as the frequencies of 
Wu modes, i.e., if we see those frequencies as functions of the stellar compactness, those are almost independent from 
the values of K. and the adopted EOS. However, in the case of wi modes, there exists the crucial different feature 
in contrast to the wu modes. That is, the dependence of the frequencies in TeVeS is obviously different from those 
expected in GR. In other words, with this different dependence of frequencies on the gravitational theory, one can 
distinguish the gravitational theory in strong-field regime by using the gravitational waves observations. In fact, the 
oscillation frequencies of axial wi modes in GR and TeVeS can be expected with high accuracy via the following 
empirical formula 



Re (ujR) = a 



f Madm \ 



(4.3) 



where (a,/3) = (2.797,4.255) in GR and (a,/3) = (2.533,3.714) in TcVcS. 

At the last, in order to discuss the possibility to make a restriction on the value of parameter JC, we show the 
dependences of frequencies of wu and wi modes on the total ADM mass in Figs. [4] and [5l These figures show that 
due to the uncertainty of EOS it is not easy to determine the exact value of K. by using the observations of axial 
Wu and wi modes. However, with the help of the information of stellar mass, combining the observations of axial 
gra vitational waves with those of frequencies of fluid oscillations, which depend strongly on the stellar average density 
[lOj . it might be possible to make a kind of constraint on the parameter K. and/or on the stellar EOS. 



V. CONCLUSION 



In order to examine how the tensor-vector-scalar (TeVeS) theory affects on the gravitational waves emitted from 
the neutron stars, we derive the perturbation equations in TeVeS with the assumption that the perturbations of 
vector field are neglected. Actually this assumption might not be physical, but with smaller value of the vector-field 
coupling constant /C suggested by the observations of neutron star, the effect due to the presence of vector-field 
perturbation on the spacetime oscillation might be too small to neglect. Thus to see the dependence of spacetime 
oscillation on the gravitational theory, we adopt this simplification as a first step. With this assumption, we find that 
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FIG. 3; For wi modes, the damping rates (left panel) and the oscillation frequencies (right panel) as functions of stellar 
compactness AIaum/R- The solid symbols correspond to the results for EOS A and the open ones to those for EOS II. 
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FIG. 4: The normalized frequencies of wn modes with two different EOSs are plotted as functions of of the total ADM mass, 
where the left panel corresponds to the damping rates and the right panel are to the oscillation frequencies. The solid lines 
denote the results with EOS A and the broken lines are results with EOS II. Additionally, the different symbols denote the 
results with different values of IC, such as the circles for K, = 0.5, the squares for K, = 1.0, and the triangles for /C = 1.5. 



the perturbations of scalar field couple with the polar perturbations of both spacetime and fluid, but they do not 
affect on the axial perturbations. Since the axial spacetime modes are known to have the same qualitative behavior as 
the polar spacetime modes (at least in the case of GR), in this article we have studied only axial gravitational waves, 
which are corresponding to the oscillations of spacetime itself. 

With two different equations of state (EOS), we calculate the complex eigenfrequencies of axial wu and w modes, 
where the real and imaginary parts of complex frequencies are corresponding to the oscillation frequencies and the 
damping rates of the emitted gravitational waves, respectively. We find that the dependences of the both frequencies 
of Wu and wi modes in TeVeS on the stellar compactness are almost independent from the parameter IC and the 
adopted EOS. The dependences of wi modes are obviously different from that expected in GR, while the dependences 
of Wn modes are similar to that in GR. Owing to these differences of dependence on the gravitational theory, one can 
distinguish the gravitational theory in the strong-field regime by using the direct observations of gravitational waves 
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FIG. 5: The normalized frequencies of wi modes with two different EOSs are plotted as functions of of the total ADM mass, 
where the left panel corresponds to the damping rates and the right panel are to the oscillation frequencies. The meaning of 
lines and symbols are same notations as Fig. U) 



emitted from the neutron stars. Additionally, we show the possibility to make a kind of constraint on the value of /C 
and/or the stellar EOS via the observations of gravitational waves. 

It should be noted that if one takes into account the effect of the vector-field perturbation, the quasinormal 
frequencies of spacetime oscillations could be changed with larger value of the coupling constant K. although the 
qualitative behavior might be similar to our results in this paper. However, if our discussion would restrict in the 
region with smaller value or /C as suggested in the neutron star observations, we could expect that the behavior of 
frequencies is almost same as that shown in this paper, i.e., we can see the obvious difference between the quasinormal 
frequencies expected in GR and in TeVeS. So still it could be possible to distinguish the gravitational theory with 
using the gravitational wave observations, although in this case it might be difficult to make a constraint on the value 
of /C. 

In this article we focus only on the axial gravitational waves, but the study of polar gravitational waves are also very 
interesting, since the gravitational waves couple directly with the perturbations of scalar field. Furthermore, since for 
simplicity we neglected the perturbations of vector field, more detailed studies are needed including the perturbation 
of this field as well as those of metric and scalar fields. Via these more detailed studies, we could obtain the additional 
information in the spectrum of emitted gravitational waves and one can provide more accurate constraints on the 
gravitational theory in the strong-field regime. On the other hand, the consideration of the effects of the stellar 
magnetic field might be of great important. In fact, recent observations of quasi-periodic oscillation in the giant flares 
of magnetar are believed to be related to the oscillations of these strongly magnetized neutron stars [13, [HI, [11| . By 
taking into account the effects of the stellar magnetic field, one might be able to set further constraints in the theory 
that describes gravity. 
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APPENDIX A: THE EFFECTS OF PERTURBATION OF VECTOR FIELD ON THE SPACETIME 

OSCILLATIONS 



If one would consider the full problem where the perturbation of every field is taken into account, the situation 
becomes more complicated. Although in this paper for simplicity we neglected the perturbation of the vector field 
such as SUi = 0, in this appendix we want to see the effect of that perturbation on the metric perturbation briefly. In 
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general, the perturbation of vector field could be given as 

'1 



-u/2 



Sgtt, SUr, 5Ue, 5U^ 



(Al) 



where the perturbed variables SlAr^ SUe, and SU^ are determined by calculating the linearized field equation for the 
vector field. With this expression of perturbation of vector field, the perturbed metric tensor in the Einstein frame is 
given by Eq. p.6p as 



/=2 m=-l 



/O SUr 5Ue 5U^ 

* 

* 
V * 



(A2) 



where Sgl^u denotes the perturbed metric tensor in Einstein frame with the assumption that 5Ui = 0, which is 
corresponding to Eq. (|3.7p . Thus if one considers the coupling of the perturbation of vector field with the perturbations 
of the other fields, the both axial and polar perturbations are connected to the vector-field perturbations. Namely, the 
axial spacetime modes should be also calculated together with the perturbation of vector field. However as mentioned 
in the main text, since the observations of neutron stars suggest that the value of the coupling constant for vector field, 
/C, should be smaller than at most K, ~ 0.8, the effect of the vector-field perturbation on the spacetime oscillations 
might be small. So as a first step we neglect the effect of the vector-field perturbation in this paper. 



APPENDIX B: THE COMPONENTS OF THE LINEARIZED EINSTEIN EQUATIONS 



In this appendix, we provide the explicit form of the various expressions used in Eqs. ()3.20p - (j3.23|) for the 
linearized Einstein equations. 



1 

+ 2 
1-K. 



1 



Ho 



v' 1 



Hi 7— V"^'/' 

k 



Pln^ - k) Je-C {^.' + C') (K - k) - \e-^k - SnGe-^'^ [p + P) 

+ ^) + Ig-C f^'^ 2u") - ^(m + rm') + (2r - 3m - rm') 
4 V J 2r^ 



2r2 



SttGc-^^ (p + P) - l6TTGPe-^f 



2r2 



2r 



Sim [Hq — H2] , 



1 ^e-^{i^'-C)hu 



(Bl) 
(B2) 



(B3) 



(B4) 

(B5) 
(B6) 



where ^ = Lp' . 



APPENDIX C: NUMERICAL TECHNIQUES 



Integrating from stellar center outward, one can get the values of $ and $' at some radial position r = Ta where 
is outside the star. With these values, we will have to match the numerical solution with the appropriate asymptotic 
boundary conditions, which is the absence of the incoming radiation. In order to find the asymptotic form of the 
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solution of Eq. I|4.ip with P = 0, we can assume a solution of the form 



$(r) 



2mr 



- 1 



(Cl) 



where rric is the mass function evaluated at radial infinity. Substituting this form of the solution into the perturbation 
equation (14. ip and keeping the terms up to the order 1 /r^ , we obtain a five-term recurrence relation for the expansion 
coefficients a„ for n > 1 



+ PnO-n + 7nan-l + 6na„-2 + <^nan-3 = 0, 

where the coefficients of the recurrence relation are given by the following formulas 

an = CQn{n + 1), 
Pn = don + cin{n - 1), 
7« = eo + di(n - 1) + C2(n - l){n - 2), 
<^n = ei + d2{n - 2) + 03(71 - 2)(n - 3), 
= 62 + dsC" - 3) + Ci{n - 3)(n - 4). 



(C2) 

(C3) 
(C4) 
(C5) 
(C6) 
(C7) 



In the above formulas, the coefficients rf^, and Ci are functions of the background quantities and have the form 

2toi 



2mc 8rni 



62 = 



C4 

C3 

C2 - 1- 

ci = -2 

Co = 1 - 



6mi 

da = 5- (2^a;TOc + 1) , 

4iujmi Grric „ , 
02 = 1 Ms, 

= 2 - 2^2 - 3^3, 
do — -2iujra - 2 + ^2 + 2^3 
24m^mia;^ Simluj 2 





ra^ ' 


6r7ic 


12mi 


fa 




6mc 
+ 


8toi 


ra 




2mc 


2toi 







— - (3mi - TO^j , 



ei 



inicUJ 



{3ml + nil) + 



262 



ep = 2mi[x)^ — 1(1 + 1) — Cl — 62, 
where mi is the coefficient in the asymptotic form of m(r), i.e., m = rric + nii/r + Oil/r^), which is given by 

kG^Ml JC 



mi 



Stt 



+ — m^. 



(C8) 
(C9) 
(CIO) 
(Cll) 
(C12) 
(C13) 

(C14) 

(C15) 
(C16) 

(C17) 

(C18) 
(C19) 

(C20) 



The first four terms of the recurrence relation (jC2ll . i.e., a_2, 0,-1, olq, and ai, are provided by the values of $ and $' 
at r = ra as 



oo 



0,1 = 



a_2 = a-i = 0, 



(C21) 
(C22) 

(C23) 
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where 

X{r) = - ij e— . (C24) 

In general, a high order recurrence relation can be reduced to a three-term recurrence relation, in which case 
a convergence criteria for the solution can be applied, and one can identify the solution describing only outgoing 
radiation (26| . In order to obtain a three-term recurrence relation, we introduce new coefficients d„, /3„, 7„, and (5„ 
as 

Oil = ai, /3i = /3i, 7i = 7i, (C25) 
"2 = Ct2, $2 = /32, 72 = 72, S2 = S2, (C26) 

and for n > 3 

otn = a„, (C27) 
/3„ = /?„-%^, (C28) 

ln=ln J , (C29) 

L^5n~ (C30) 



Then the original five-term recurrence relation \C2\ becomes a four-term recurrence relation for n > 1 

d„a„+i -I- /3„a„ -I- 7„a„_i -I- 5„a„_2 = 0. (C31) 
Furthermore, defining another set of coefficients a„, /3„, and 7„ as 



and for n > 2 



(C32) 

(C33) 
(C34) 

(C35) 

the four-term recurrence relation (IC31|) can be reduced to a three-term recurrence relation for 71 > 1 

Q;„a„4.i -I- ^„a„ -f 7„a„_i = 0. (C36) 

Using this three-term recurrence relation, the boundary condition at radial infinity can be expressed as a continued 
fraction relation between q:„, /?„, and 7„ as 



h = 


Pi, 


7l = 


OLn = 


an, 




Pn- 




l5n 


In 


-1 


% - 


k- 


l5n 


In 


~1 



ao pi- p2- P 
Since this relation can be rewritten as 



ai -71 ai72 ^273 _ ^ ^^^^^^ 



„ 3 «07l «172 0:273 _ rf X fnoQ\ 

O^Po-j-j-j-- = fi.), (C38) 

where Pq = ai/ag and ao = —1, hence, one can determine the eigenfrequnecy w of a quasinormal mode by solving 
the equation f{uj) = 0. 
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